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Abstract 

We consider the issue of an accurate description of the evolution of the non- 
singlet structure function moments M„(Q) near heavy quark threshold. To this aim 
we propose a simple modification of the standard massless MS scheme approach to 
the next-to-leading QCD analysis of DIS data. We apply it to the processing of the 
modern CCFR data for xF^ structure function and extract the value of 

UsiM^) K, 0.108 ±0.004 



We check also the consistency of light gluino hypothesis with CCFR data. 



1 Introduction 



An important means of verification of the validity of pQCD (that is, of perturbative QCD 
"improved" by the RG summation) is an analysis of deep inelastic scattering (DIS) data. 
To interpret these data within pQCD, one should pay credit to a number of subtle physical 
effects: contributions of high twists, nuclear effects, high-order (three-loop) corrections and 
the influence of thresholds of heavy particles. All the introduced corrections are roughly of 
the same order of magnitude. 

This paper is devoted to the problem of influence of the thresholds of heavy quark (HQ) 
on the pQCD analysis of DIS data that includes, in particular, the evolution of the strong 
coupling constant as{Q). Recent estimates performed in [|I], ^, ^ have revealed a significant 
role of threshold effects in the cXs{Q) evolution when the DIS data lie close to the position 
of "Euclidean-reflected" threshold of heavy particles. The corresponding corrections to 
as{Mz) can reach several per cent, i.e., they are of the order of the three-loop and 
nuclear effects on as{Mz) ■ 

A common algorithm for the renormalization-group (RG) resummation is based upon 
beta-function /9(as) and anomalous dimensions 7(0:5) calculation and the RG differential 
equations integration performed within the MS renormalization scheme. However, the 
widespread massless MS scheme fails to describe the data near thresholds of heavy particles 
- 6, c quarks and, maybe, light superpartners 0|. 

An appropriate procedure for the inclusion of threshold effects into the Q^-dependence 
of OLsiQ) in the framework of the massless MS scheme was proposed more than 10 years ago 
1^, : transition from the region with a given number of flavors / described by massless 
as{Q', f) to the next one with / + 1 ("transition across the Mj+i threshold") is realized 
here with the use of the so-called "matching relation" for as{Q) 0- The latter may be 
considered as the continuity condition for ds{Q) on (every) HQ mass 

as{Q = My+i; /) = a,(Q = M/+i; / + 1) (1) 

that provides an accurate as(Q)-evolution description for Q values not close to the thresh- 
old region. The condition (|ip is used up to the three-loop level; the other version of the 
matching can be found in [Q. 

One needs also one more element, the matching procedure for the evolution of the struc- 
ture function moment M„(Q, m). The corresponding expressions for anomalous dimensions 
7(j)(n; /) are well known in the MS scheme for a fixed / value (see, e.g., [^) but until now 
there is no recipe for obtaining a continuous interpolation across the HQ threshold for the 
moment evolution. 

In this paper, we are going to focus just on this aspect of the problem: how does the 
HQ threshold influence the evolution of the DIS structure function? We will examine only 
non-singlet processes of DIS so as to pass over the dehcate problem of modification of the 
operator product expansion in DIS through introducing a new scale, the mass of a heavy 
parton (for discussion, see 0). 

To solve the problem we propose a rather simple modification of the massless MS scheme 
to take into account thresholds in analyzing the moments of the DIS non-singlet structure 
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function at the two-loop level. 

To simplify the exposition, we shall take advantage of the explicit analytic mass- 
dependent RG-solution derived in and ||T0| that is expressed directly in terms of 
asiQjfn) and Mn{Q,m) perturbation expansion coefficients. This allows us to avoid the 
use of RG-generators, that are [3 and 7-functions. In the next section, we present smooth 
analytic expressions for the evolution of dsiQ, ^) and M„((5, m) at the 2-loop level based 
on this mass-dependent RG formalism. We shall omit all theoretical and technical details 
(they can be found in refs. 0, |11|) and write only final results. In Sect. 3, we introduce 
the "spline-approximation" to describe the two-loop level continuous moment evolution, and 
present there another proof of the matching condition ([^). In Sect. 4, we describe briefiy 
a method of analysis of the DIS data. On its base we carry out the fit of fresh CCFR 
Collab. data, extract the parameter as and estimate the contribution of threshold effects. 
We discuss the consistency of MSSM light gluino existence with CCFR data by using the 
spline-type evolution of M„((5,m) in Sect. 5. 

Throughout the paper we use the notation: a = as/4:7i, (a = as/47r); indices in brack- 
ets stand for the loop number, e.g., = Pe-i] instead of the structure function moments 
Mn{Q,m), we consider only its "evolution part" A^„((5,m) (i.e., moments of the distribu- 
tion function) 

Mn{Q,m) = Cn{a,Q,m)-Mn{Q,m), (2) 
where Cn are moments of the coefficient function of a certain DIS process (see, e.g., 



Ref. m) 



2 Mass-dependent RG solutions 



In the massless case, the moment two-loop evolution is described by the expression 
with the numerical coefficients 



M^^\Q) = Mn{^^) i-^A exp{[a-a(^)(Q)]/„} (3) 



Po Po 
/3o = = 11 - Ml)/ ' Ml) = 2/3 ; A = A2)(/) = 102 - A/3(2)/, A/?(2) = 38/3 . 

In the mass- dependent case, one should use instead of Eq.(||), a bit more complicated 
expression p of the same structure 

Dn{Q,...) 

with the functional coefficients Dn, Fn 



M^:\Q) = MM " ) exp { [a - a(2)(Q,m)] F„(Q, ...)} (5) 



Dn{Q,m,^) = .^^V/f'*^^ ^ ; (6) 
A(i)(Q,m,/i) 
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Fn{Q,m,fi) 



(g, m, /i)r(2) (n, Q) - A(2) (Q, m, /i)r(i) (n, Q) 



[A(i)(Q,m,/i)]2 
and the two-loop running coupling a taken in the form 



(7) 



a(2)(Q, rn; a) a ' 
In the non-singlet case of DIS, 



v4(i)(Q,m,/x) + 



A(i)(Q,m,/i) 



ln[l + aA(^i){Q,m, jj,)] 



(8) 



1(1) (n, m, Q) = 1(1) {n,l) = 7(1) (n)/, / = In (qVK 

and the HQ-mass-dependent y4(^), r(£=2)('^, Q) appearing in Exp. (H-H) are just perturba- 
tion expansion coefficients: 



a{Q,m,n;a)pert = a- a A(^i){Q,m,n) + a-^ < A^i){Q,...) - A(2)((5, m, /i) ^ +...; 



MniQo 



1 + ar(i)(?2, /) + 



pert 



r(i)(n,/) (r(i)(n,/) - A(i)(Q,m,/i)) 
+ a < ^ 7: ^ + r(2)(n,Q,m,/i) ^ + 



(9) 



satisfying the normalization condition - v4(^)(Q = /x) = r(£)(n, Q = yu) = 0. 

These coefficients consist of the usual massless part (for / = 3) and HQ-mass dependent 

contributions, e.g., 



A(i)(g,m,/i) 



1 



1 



/?(i)(3)/-A/?(i)^ 



mi 



(10) 



with summation over HQ's: h > A. Here Ji is the one-loop fermion mass- dependent contri- 
bution, like the polarization operator or the three-gluon vertex loop |T^ subtracted 
at = and C being some subtraction scheme-dependent constant. 

"Massive" RG solutions (^ and (§) possess several remarkable properties: 

• they are built up only of "perturbative bricks", i.e., loop-expansion coefficients ^(£), 
r(£)(?T,) (taken just in the form they appear in the perturbative input) and "contain 
no memory" about the intermediate RG entities such as j3 and 7 functions; 

• in the massless case with pure logarithmic coefficients, A(^k) = (3{k)h '^(k){n) = 7(A;)(n)/, 
they precisely correspond to the usual massless expressions, like Eq.(||); 

• being used in QCD, they smoothly interpolate across heavy quark threshold between 
massless solutions with different flavors numbers. 
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3 Smooth schemes and MS massless schemes 



3.1 Smooth mass-dependent scheme 

We have above considered general formulae to describe the A^n(<5)-evolution including 
the threshold effects. It is clear that the mass-dependent MOM schemes automatically 
provide the most natural smooth description of thresholds. To use them in the framework 
of leading order, one needs a mass- dependent expression for Ji presented of Appendix A 
(for two different schemes). So, to perform the one-loop evolution analysis of moments, one 
should substitute Eq. ( p3D or (pTl) in Appendix A into Eq. (|T0|) and then into Eq. (Q) and 
Eq. d^), and use the approximation 



We shall perform the fit of CCFR data following this formula in the next section. 

However, the MOM scheme meets tremendous calculational difficulties in the next-to- 
leading order of pQCD. Moreover, each of the expansion coefficients y4(j), r(j) (see Sec. 2) 
becomes gauge-dependent at the two-loop level, which is not convenient. These difficul- 
ties are absent in the widespread MS scheme. One can go far in loop calculations here 
(see |12|), but the scheme is not sensitive to the thresholds at all. Below we suggest 
a practical compromise between these different possibilities - the "spline" scheme. This 
scheme possesses both the sensitivity to thresholds and simplicity of the MS procedure. 
Nevertheless, the MS scheme looks like a conventional standard for all DIS calculations 
now. Therefore, one should recalculate the results obtained in other schemes to the MS 
scheme at an appropriate number /. We do not need recalculation for a^^'^^^Mz) because 
the spline and MS (at / = 5) schemes evidently coincide at /x = Mz by construction, i.e., 
afSspi^Mz) = af^{Mz-J = 5). 



3.2 MS — vulgate scheme 

Usually, to obtain the evolution law, one calculates numerical expansion coefficients of 
generators P{as), 7„(as)... in MS scheme and solves the massless RG equations. In the 
solution, with all integration constants being omitted, one arrives at the final procedure 
which we shall name the MS -vulgate scheme. 

The first recipe to include the threshold mass Mj into the framework of the MS evolution 
was formulated in Ref. as the "matching condition", Eq. (|1]), for the coupling constant. 
Now all measurements on a low scale Q are usually interpreted in terms of the as{Mz) 
- RG solution in a certain scheme, with an appropriate matching of different numbers of 
active flavors which evolve from the scale Q to Mz- The matching condition (|l|) leads to a 
simple rule of including next "active flavors" h into the evolution law 




(11) 



^fff(0=%(3)^-A/3»-W - 



(Q,...) = /5(,)(3)/-A/3(,)r; 



(12) 
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i* = i: [oiQ' - in(g/M,(,))2 - (g ^ ^,) 

h 

Nevertheless, the threshold value M/i(j) does not follow from this procedure and is left 
uncertain. 

Note, the spline-type (in terms of the /-variable) expression (|T2p has an evident anal- 
ogy with the approximation for the mass- dependent MOM scheme formulae for 74(j)(g,yu) 
( r(j)(n,g,yu)) with the structure y4(i)(g,/i) ~ (/(i)(g,m) - m)) , see, e.g., Eq. 
The approximation being discussed needs an asymptotic form of the mass- dependent cal- 
culation for the elements I{i){z = Q'^/m?), i.e., only logarithmic and constant terms 
I[i){z) — > ln(2;) — Cj (see, e.g., Exp. (p^ and (^) in Appendix A). Based on this form 
one can construct a simple "pure log" ansatz for /(j)(g,m): 

/(,)(g,m) ^ jg^^^^P'(g,M,(,)) = e{Q^ - MI^^) in (gVM'w) ; = m,exp(Q/2). 

This ansatz roughly imitates the "decoupling" property of I(^i){Q,m) at g < M^-) and 
provides its asymptotic form at g > It leads to the approximation for y4(j)(g,/i): 

A(,)(g,/i)^ 4f^^^^'(/) = 

/?«(3)/ -A/5« ^ [^(g^ - M,2(,)) In (gVM.%) - (g - /i)] , (13) 

where the threshold position M/j(j) is determined by the scheme dependent constant C(j). 

A certain value of the threshold M/i(j) in Exp. (12) and another proof of the matching 
condition (|I|) for the MS scheme can be obtained by using, e.g., the "three-step procedure" 
introduced in p|. Let us review it briefly. 

Well below the threshold, for yU <^ M, one usually uses some effective MS scheme, 
say MSi, that does not take mass of a particle into account. Above the threshold, a new 
particle cannot be ignored, but when Q ^ M, it can approximately be treated as massless 
within some other MS2 scheme. How should the couplings ai(/i) and 02 (g) in these two 
MS schemes be related? The answer can be obtained by the three step algorithm: 

(i) recalculating from the MSi to MOM scheme at g = ;U < M to get a*^'^^^(/i); 

(ii) performing the RG evolution of 

^MOM up to g = g > M in the MOM scheme; 

(iii) recalculating to the MS2 scheme, including the mass contribution, at g = g. 

The final result of these successive steps leads to the approximate (due to power corrections) 
equality Mh ~ rrih [@] for the threshold at the two-loop level. We obtain just the same 
result as usually used for the matching condition mentioned above with Mh = rrih. 

Consequently, proceeding in this way, one must modify the perturbative expansion 
coefficients for ,i.e., r(^{^{n, Q) in (^), in the same manner as the expansion coefficients 

of the coupling constant y4(j) = A^^'^^K For this aim we recall the structure for 

7(2) ('^; /) in the framework of the MS scheme (for details see, e.g., |]T^) 

7(2) (n; h + 3)= 7(2) (^; 3) - A/?i ■ /i ■ A7(2)(n); (14) 
A7(2)(n) = y |l0^i(n) - 6^2(71) - ^ - ; Sk{n) = ^, 
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where the first term 7(2) (n; 3) in the r.h.s. of Eq. (JT^) consists of the usual massless 
part and the parameter h numbers here heavy flavors. It is known (see e.g. |16|) that 
7(2) (^;/) contains the terms generated by the evolution of the coupling constant. These 
terms naturally appear in the calculation of two-loop diagrams for 7(2) ('t-;/), they are 
proportional to the coefficient (see the second term in (0)). Therefore in the MS - 
expression for r(2) there appears a term proportional to the one-loop coefficient ^(i): 

1(2) (n, Q;h + 3) = 1(2) (n, Q; 3) - A(3,{h ■ I) ■ A7(2)(n); (15) 

1(2) (n,Q; 3) = 7(2) (n, 3)/. 

The /i-dependent part of the A(i)-term is singled out of Exp.(|l5]) in the form of A(3i{h ■ I). 

To obtain the continuous coefficient r^^^''^^(n, Q, h+3), one should substitute ^4(1) A^^'''"'^'', 
i.e., {h ■ I) I* into the second term of the r.h.s. of Eq. ([T5|), according to the recipe 

(0): 

1(2) (n, Q;h + 3)^ r^'^^n, Q-h + 3) = 1(2) (n, Q; 3) - Af3^l* ■ A7(2)(n). (16) 



Now we can get the complete evolution law by substituting (12) and (|T^) into formulae (H) 
and (^, (0) and then into the general formula (||): 



/ \ D„{Q,M) 



exp 



Recent CCFR CoUab. experimental data on DIS are processed by this method (taking also 
account of the one-loop coefficient function) in the next section. 



4 The QCD fit of the xFs CCFR data 
4.1 Method of QCD Analysis 



In this section, we present the QCD analysis of the CCFR data [|T8[. They are the most 
precise data on the structure function xF^^XyQ^) . This structure function is pure non- 
singlet and the results of analysis are independent of the assumption on the shape of gluons. 
To analyze the data, the method of reconstruction of the structure functions /^from their 
Mellin moments is used |1^ . This method is based on the Jacobi - polynomial expansion 
of the structure functions. 

m I 



Following the method 



we can write the structure function xF'? in the form: 



N„ 



{x,Q' 



X 



XI 



n=0 



X 



3=0 
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where G"^(x) is a set of Jacobi polynomials and c"(a,/3) are coefficients of their power 
expantions: 

n 

e:'^{x) = Y.4'\a,(3)x^. (19) 

j=0 

The quantities A'maa', « and (3 have to be chosen so as to achieve the fastest convergence 
of the series in the r.h.s. of Eq.(^) and to reconstruct xF^lxjQ"^) with the accuracy 
required. Following the results of |jT9| we have fixed the parameters - a = 0.12 , f3 = 2.0 
and Nmax = 12 . These numbers guarantee an accuracy better than 10"'^ . 

Finally, we have to parameterize the structure function xFs{x, Q'^) at some fixed value 
of = Qq ■ We choose xF^i^x, Q"^) in a little bit more general form as compared to pT 



where the same data have been analyzed within QCD in terms of A^^^ without thresholds 
effects: 

xF3(x,Q^) = Ax^(l-x)^ (l + 7x). (20) 

Here A, B, C, 7 and = ds{Qo) are free parameters to be determined by the fit. 

To avoid the influence of higher-twist effects, we have used only the experimental points 
in the plane (x, Q^) with 5 < < 501 {GeV/cy and 0.015 <x< 0.65 . The effect of 
target-mass corrections in xF^ is taken into account to order M'^/Q'^ ^ 



4.2 Results of Fit and Discussion 

Here we present the results of processing the CCFR data obtained in the framework of two 
different approaches: 

First, we have used the massless /-fixed MS - scheme approach based on the two- loop 
evolution formula @. The corresponding results for ao{f) are collected in the up-part of 
Table 1; 

the second approach is based on mass- dependent evolution Eq.(H), for this formula we 
adapt the "spline" approximation (|1^. These results for af^^ are collected in the down-part 
of Table 1. 

Both the parts of Table 1 include the results of the LO and NLO fit; they are completed 
with the results of evolution of q;o(/) and af'' to the point Q = Mz with appropriate 
matching (P of different numbers of active flavors (the last column - as{Mz)', see, for 
comparison with other estimations, recent review p4|)- Repeating the fit procedure for 
different values of Ql = 5, 50, 500 GeV^ we have obtained the experimental dependence of 
ao{f) on the momentum transfer. In all fits, only statistical errors are taken into account. 
Here we make some comments on the fit results. 

• A significant decrease of ;i^2(A^iO) comparison with ■x^^^'^'> for all variants of the fit 
in Table 1 demonstrates that second loop effects are important for the kinematical 
region under consideration. 

• To demonstrate the sensitivity of ao{f) on the /, the results of the fit both for LO 
and NLO are shown for different / in the up-part of Table 1. The largest difference 
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between the values of ao(/) for different / is about 5% in the case of Ql close to kine- 
matical boundaries Qq = 5 and 500 GeV"^. This difference reduces to 2% variations 
for as{Mz] /). There are opposite relations for ao(/) for these two points: 
ao(3) > a;o(4) > ao(5) for Ql = h GeV"^ and a;o(3) < a;o(4) < a^ih) for = 
500 GeV^. Note, the effects mentioned above can be described by estimating Aa^ = 

l( dfMnia,f,Qexp) \ 
"''^^V daMn{aJ,Qexp)) 

Here, the brackets ((...)) denote the average over experimental values of Ql^p- At the 
one-loop level, this expression leads to the simple estimate 

^~-A/5(i)a.((ln(^))), 

in qualitative agreement with the results in Table 1. 

• The final results for as{Mz) depend on the Ql choice. For the /-fixed scheme this 
dependence amounts to 2% (two thousandth of the absolute value) and for the spline 
scheme, it is within only 1% (one thousandth) of the value of as{Mz)- Of course, 
the value of af\Qo) lies between the corresponding values of (y^^{Qo',f = 4) and 
a^^{Qo; f = 5). This suppression of the residual dependence on Ql for the spline 
scheme results gives an additional "phenomenological" hint for preferring this scheme 
over the MS /-fixed version. 

• The results of the fit are rather stable to the mass variations. The 10% change of 
and Mfe yields less than 0.5% change for a^. 

Comparing the up- and down-parts of Table 1, we arrive at two main conclusions: 

1. The spline scheme is more preferable than the traditional massless scheme, to process 
the experimental data involving thresholds, the values of al^\Mz,Qo) are focused 
tightly. The average value a^^^Mz) = 0.108 is considerably smaller than the average 
LEP result and about 1 s.d. lower than that the CCFR result for as{Mz) obtained 
recently ||2^ from the Q^-evolution of xFs and F2 for > IbCeV^ {as{Mz) = 
0.111 ± .004). 

2. The threshold effects reveal themselves as approximately a 1 % correction to the value 
of as{Mz)- The final results for the average values of as{Mz) in the fit presented in 
Table 1 look like: 

af{Mz) = 0.108 ± 0.001 (t/ieor) ±0.003 (stat) 
ap(Mz;/ = 4) = 0.107 ±0.002 (theor) ±0.003 (stat). 

Theoretical errors presented here include the uncertainties due to Jacobi polynomial tech- 
nique reconstruction and Qp-deviation of the as{Mz) value in the fit. 
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As it was emphasized in Sect. 3.1, the most natural way to include the thresholds effects 
into the analysis of data is to use the MOM scheme in formula @. We present the results 



of processing the same experimental data by the one-loop evolution formulae (|TT]) in Table 
2. Note that one cannot directly compare the results obtained by different schemes, the 
spline fit or the MS fit in Table 1 and the MOM scheme fit in Table 2 (see Sect. 3.1). 
Indeed, we extract, generally speaking, different kinds of the coupling constants, 
^MOM ^ To connect these quantities, we need relations between the MS and our version of 
MOM-scheme. Here we write the relation at the one-loop level on scale ii 

^ = j;Mc^ + Oa)M + --- (21) 

The coefficients depend, in general, on the gauge parameters and on the choice of a 

particular MOM scheme. The expression for was represented, e.g. in [|l^, with the 

5 

constant Di = - ■ /3(i)(3) — (11/12)C^ in the Landau gauge: 
3 

Q(0(/iV^L ••) = ^ (em' ■ (W/rnl) - ln(/iVm^j) + a) • (22) 

So, the quantity a^^{Mz) in the second rows of Table 2 was obtained by successive sub- 
stitutions Ji into Eq.(P^ and than into Eq.(pTD. 



5 CCFR data and the light gluino window 

In Subsec. 5.2, we have obtained a comparatively small effect, one-two thousandth to 
as{Mz) for the threshold contribution at the fit. The order of the effect is determined by the 
reason that only one 6-quark threshold "works" really in the region in question. If Nature 
would provide few thresholds in the experimental interval 5 GeV^ < < 501 GeV^, 
then they combined influence in a fit becomes significant and the preference of the spline 
scheme should look evident. To demonstrate this here we took an attempt to reconcile the 
existence of light MSSM gluino (g) and the CCFR data. The possibility of the light gluino 
existence {rrig is of the order rrib) was intensively discussed few years ago in the context 
of the discrepancy between low energy as values and the LEP data at the Mz peak. This 
discrepancy has a chance to be resolved by including the light gluino The Majorana 

gluino leads to large effects in the evolution - A/S^^s^ = 2 , = 48 in Eq. (^ and slows 

the running both the coupling constant, and the moments AiniQ)- This reinforcement of 
the contribution of a new g threshold must influence the fit parameters. It is clear that the 
standard MS -scheme at a fixed / everywhere is not adequate to the situation. We have 
performed the fit of CCFR data for different values of the gluino mass rrig, the results of 
the fit for at rrig = 3, 4 GeV and rrig = 10 GeV are presented in Table 3. 

Rather a strong dependence of a*^'(^^(M^, Qq) value on Ql (three thousandth), as well as 
a slight growth of x^, do not provide confidence that the data are consistent with the gluino 
with mass rrig < 4 GeV (see up-part of Table 3). The "non-renormgroup" dependence of 
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af'-^'^^MzyQo) on Ql reduces less than 1.5 thousandth for nig > 10 GeV. This variation 
is not too far from the variation of in Table 2, so the gluino with these masses may be 
considered as accessible for the data. This leads to the estimation for as{Mz)'- 

at m-g = 10 GeV, af^^\Mz) = 0.119 ± 0.002 {theor) ± 0.004 {stat), 



that is close to the LEP data P3. Nevertheless, if one considers clS cl fit parameter and 



"releases" it, the best x will be reached at nig beyond the fitting region, i. e. 



m-g > 22 GeV. 



6 Conclusion 

We have devised here a new approach for describing the two-loop evolution of the moment 
AiniQ"^) of structure function of lepton-nucleon DIS involving the threshold effects of heavy 
particles. The approach employs an analytic quasi-exact two-loop RG solution 
within the framework of mass- dependent Bogoliubov Remormalization Group (see Refs. p^ , 
p7[]), and on the results of papers [0, flj. 

We adapt here the spline approximation ( p!7|) to the above mentioned RG solution 
and obtain as a result the simple modification of the formulae of the standard massless 
MS evolution. For the particular case of a coupling constant evolution this approximation 
effectively leads to the same result as the "matching condition" . Finally, this recipe provides 
a more realistic continuous description for the A1„((5^)-evolution and looks rather simple 
from a practical point of view. We performed the processing of the modern CCFR data to 
extract the value of as{Q) by three different ways: 

(z) the traditional MS-scheme at the fixed numbers of fiavors /; 

(m) the spline scheme with break point at mass M = m; 

{iii) the MOM scheme (in the leading approximation); 
and compare results of fits. 

The results for the MS spline - scheme processing of the data are the most adequate to 
the situation both for the physical and practical points of view. The threshold contribution 
to the value of as{Mz) consists of about 1%; the extracted value of af\Mz) is equal to 
0.108 ±0.004. We examine the possibility to reconcile the CCFR data and the MSSM light 
gluino. It is possible for gluinos with mass nig > 10 GeV, but it seems that similar gluinos 
are less probable due to other constraints (see P5|). 
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Note Added in Proof 

After this work was completed and submitted to the hep-ph database [a], the CCFR col- 
laboration announced in Ref. [b] that the corrections of the energy spectrum calibration of 
the neutrino beam in the CCFR experiment change the "old" CCFR data of Ref. [18] and 
lead to the essentially larger values of as{Mz) '■ 

xFs : a,{Mz) = 0.118 ± 0.0025(stot) ± 0.0055{syst) ± 0.004(t/ieory) 

extracted at the NLO of perturbative QCD using the integro-differential approach based 
on the CLAP equation. It would be better to verify this result by other methods. We 
believe that our qualitative result on the order of the "threshold correction" to the value 
of ds{Mz) could not be changed by any modification of the experimental data and would 
be the same ~ +10~^. 

Additional References 

a. D.V. Shirkov, A.V. Sidorov, S.V. Mikhailov, |hep-ph/ 9607472 

b. CCFR/NuTeV CoUab., D.Harris, talk at the XXVIII Int. Conf on HEP, Warsaw, July 
1996. 
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A Appendix 

Here we list the mass- dependent one-loop expressions for Ii{z) and the corresponding f3 
functions in the MOM scheme. First, we write the well-known exact expression for the 
fermion polarization operator ( z = /mf ) 

Il^{z) = 2^1 + A/z (l - ^) In (v^l + z/4 + y^) + ^ - ^ ; (23) 

If{z ^oo) = \nz-{c, = ^) + 0(1/ z) ■ Il{z - 0) = ^ + 0{z^) . (24) 
and the {3 function 

/3(i) = ^c^-|t,x:(i-/^.W); 

•J i=l 

where hi{z) = h1{z) = -^- ^3^, ,^ In {^l + z/A + (25) 



1 + 2/4 

The three-gluon contribution Iq{z) can be expressed in terms of dilogarithms Li2{z), the 
final formula for Iq{z) is too cumbersome and we do not demonstrate it here. As it was 
predicted in |T^, the result for the corresponding three-gluon contribution to the /?- 
function, hi{z) = hf{z) -function may be expressed in terms of dilogarithms, as well. 



1 



-6<^ — , In U/1 + z/A + J z/A + 

1 v/i(i + 2/3) yi^T/i / V / / 

+ --^Im [Li2{z*^) - Li2{z3) + Li2{zl) - Li2{zi)] 



I 



V ^ v^/ \2^l + z/A + i^) ' V V^) \2^l + z/A-i^^ 

The function hl{z) is described by the rational approximation |T5| : 

This approximation works well (better than 1% of accuracy) when 2; > 1, but when z k, 
10^^ — 10^^, the accuracy is about 10%. One can restore the corresponding approximate 
expression for Iq{z) by the elementary integration of h^{z) 

Ifiz) = A ■ ln(l + z/zi) + (1 - A) ■ ln(l + z/z2) ; (27) 
A = 21/10 ; zi = 20/3 ; Z2 = 5/2 ; 
If{z^ 00) = In z -{cg ^2.98) + 0(1/ z) ; (28) 



12 



References 

[1] Yu.L. Dokshitzer and D.V.Shirkov, Z.Phys.C 67 (1995) 449. 
[2] D.V. Shirkov and S.V. Mikhailov, Z.Phys.C 63 (1994) 463. 

[3] J. Chyla Phys. Lett B351 (1995) 325; A.L.Kataev et al., INP-0919/96, Moscow, 
1996; JINR E2-96-91, Dubna; hep-ph 9605367 

[4] A.V. Sidorov, M.V. Tokarev Phys. Lett B358 (1995) 353 . 

[5] W.Bernreuther and W.Wetzel, Nucl. Phys. B197 (1982) 228. 

[6] W. Marciano, Phys. Rev. D29 (1984) 580. 

[7] A. Gonzales-Arroyo, C. Lopez, F.J. Yndurain Nucl. Phys. B153 (1979) 161. 

[8] B.J. Edwards and T.D. Gottschalk, Nucl. Phys. B196 (1982) 328. 

[9] D.V. Shirkov, Theor. Mat Phys. (USA) 49 (1981) 1039. 
[10] D.V. Shirkov, Nucl. Phys. B371 (1992) 267 . 
[11] D.V. Shirkov, JINR Preprint E2-95-341. 

[12] S.A. Larin, T. van Ritbcrgcn and J. A.M. Vcrmaseren, Nucl. Phys. B427 (1994) 41. 

[13] F.J. Yndurain, Sec. 23 in "Quantum Chromodynamics. An Introduction to the The- 
ory of Qarks and Gluons"( Springer, 1983). 

[14] H. Georgi and H.D. Politzer, Phys. Rev. D14 (1976) 1829. 

[15] O.Nachtmann, and W. Wetzel, Nucl. Phys. B146 (1978) 273. 

[16] S.V. Mikhailov and A.V. Radyushkin Nucl. Phys. B254 (1985) 89. 

[17] T.Yoshino, K.Hagiwara, Z.Phys.C 24. (1984) 185. 

[18] CCFR CoUab., S.R.Mishra et al., Nevis Preprint N 1459 (1992); CCFR CoUab., 
W.C.Leung et al., Phys. Lett. B317, 655 (1993); CCFR CoUab., P.Z.Quintas et al., 
Phys. Rev. Lett 71, 1307 (1993). 

[19] V.G.Krivokhizhin et al, Z. Phys. C36, 51 (1987); C48, 347 (1990). 

[20] G.Parisi, N.Sourlas, Nucl. Phys. B151, 421 (1979); I.S.Barker, C.B.Langensiepen 
and G.Shaw, Nucl. Phys. B186, 61 (1981). 

[21] A.L.Kataev and A.V.Sidorov, Phys. Lett B331, 179 (1994). 

[22] O. Nachtmann, Nucl. Phys. B 63 (1973) 237; H. Georgi, H.D. Politzer, Phys. Rev. 
D 14 (1976) 1829; S. Wandzura, Nucl. Phys. B 122 (1977) 412. 

[23] CCFR-NuTEV CoUab., D. Harris et al, FERMILAB-Conf-95-114, March, 1995. 

[24] S. Bethke, Nucl. Phys. (Proc. Suppl.) 39 B, C (1995) 198. 

[25] J. Ellis, D.V. Nanopoulos, D.A. Ross, Phys. Lett B305 (1993) 375; U.Amaldi, W. 

de Boer and H.Furstcnau, Phys. Lett B260 (1991) 447; I. Antoniadis, J. EUis, D.V. 
Nanopoulos, Phys. Lett B262 (1991) 109; M. Jczabck and J. Kuhn, Phys.Lett.B 301 
(1993) 121; J. Ellis, D.V. Nanopoulos, D.A. Ross, Phys. Lett B305 (1993) 375. 



13 



[26] N.N. Bogoliubov and D.V. Shirkov, Chapter Renormalization Group in monograph 
"Introduction to the Theory of Quantized Fields", Wiley- Interscience, N.Y., 1959, 
1980. 

[27] D.V. Shirkov, Bogohubov Renormalization Group (second edition) Preprint JINR 
E2-96-15, Dubna 



14 



Quark 
content 




1-loop 






2-loops 












ao 


asiMz) 


MS - scheme, / - fixed, Mg = 


Ql = ^GeV^ 


uds 


104.4 


0.2629±0.019 


0.1097 


84.7 


0.2327± 0.014 


0.1085 


udsc 


103.5 


0.2557±0.019 


0.1133 


85.3 


0.2282 ±0.013 


0.1076 


udscb 


102.7 


0.2498±0.016 


0.1173 


86.1 


0.2238±0.013 


0.1066 


Ql = 50Ge\/2 




uds 


103.2 


0.1850±0.009 


0.1103 


84.0 


0.1648±0.007 


0.1069 


udsc 


102.6 


0.1856±0.009 


0.1139 


84.7 


0.1668±0.007 


0.1077 


udscb 


102.0 


0.1860±0.009 


0.1177 


85.4 


0.1688±0.008 


0.1085 


Ql = 500GeV^ 


uds 


104.2 


0.1421±0.0053 


0.1105 


87.3 


0.1275±0.0035 


0.1038 


udsc 


103.6 


0.1448±0.0056 


0.1140 


87.9 


0.1311±0.0043 


0.1061 


udscb 


103.1 


0.1474±0.0055 


0.1177 


88.5 


0.1348±0.0046 


0.1085 




Spline, Mc = 1.3 GeV, Mb = 5 GeV 


af{Mz) 


Ql = 5GeV^ 


uds+(c+b) 


102.2 


0.2526±0.017 


0.1167 


85.8 


0.2262±0.012 


0.1072 


Ql = 50GeV^ 




uds+(c+b) 


102.1 


0.1856±0.009 


0.1175 


85.1 


0.1682±0.007 


0.1083 


Ql = 500GeV^ 


uds+(c+b) 103.4 


0.1473±0.0058 


0.1176 


88.2 


0.1344±0.0046 


0.1082 



Table 1 The results of the LO and NLO QCD fit in the MS-scheme of the CCFR 
r5[ xFs structure function data in a wide kinematical region: 0.015 < x < 0.65 and 



5GeV^ < Q^ < bOlGeV"^ {Nexp.p. = 81). The value of the coupling constant is determined 
for different numbers of the flavor (in up-part of Table) and for few values of the momentum 
transfer ao = a{Ql = 5, 50, 500 GeV"^). We present in the first column the flavor 
content involved in QCD-evolution using in fit. The results of fit with the matching at the 
thresholds corresponding to rric = l.SGeV and nib = 5GeV are presented in the down-part 
of the table. The value of as{Mz) is calculated with a statistical error about ±0.003. 
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Quark content 








MOM-schcmc, Ql = 3GeV 


2 


uds+c+b+t 


103.4 


0.281 + 0.018 


0.115 + 0.0031 


MOM^ 


MS Ql^SGeV 






(f=4) 0.282 + 0.022 


0.106 


MOM-scheme, Ql = lOGeV^ 


uds+c+b+t 


103.1 


0.230 + 0.014 


0.116 + 0.0037 


MOM ^ 


MS Ql^lOGeV^ 




(f=4) 0.188 + 0.014 


0.106 


MOM-scheme, Ql = bOOGeV^ 


uds+c+b+t 


104.3 


0.1455 + 0.0057 


0.116 + 0.0035 


MOM^MS Ql = 500GeV 


2 


(f=5) 0.130 + 0.0052 


0.106 



Table 2 In the first rows of the Table 2 the results of LO MOM-scheme fit for different 
values of Ql are represented. The results of its transformation from cc^'^^ to oi^^{f — 4) 
or a^^{f — 5) are located in the nexts rows. 
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CCFR data and the light gluino window 



Particle 




2-loops 




content 








Spline, Mc = 


L.3 GeV, Mb = 5 GeV, M^ = 3 GeV 




uds+(c+b+5f) 


88.4 


0.2162±0.012 


0.1230 




Ql 


= bOGeV^ 




uds+(c+b+^) 


87.2 


0.1744±0.008 


0.1251 


Spline, Mc = 1.3 GeV, Mb = 5 GeV, = 4 GeV 




Ql 


= 5GeV^ 




uds+(c+b+^) 


87.8 


0.2188±0.012 


0.1221 




Ql 


= 50GeV^ 




uds+(c+b+^) 


86.7 


0.1739±0.008 


0.1248 




Spline, Me = 1.3 GeV, M^ = 5 GeV, Mg = 10 GeV 




Q't 


= 5GeV^ 




uds+(c+b+^) 


85.8 


0.2254±0.013 


0.1185 




Ql 


= 506-6^^2 




uds+(c+b+^) 


85.7 


0.1688±0.007 


0.1201 




Ql 


= 500^6^2 




uds+(c+b+^) 


89.1 


0.1392±0.005 


0.1186 



Table 3 The results of fit with the matching at the thresholds corresponding to rric = 
1.3GeV, iTLb = 5GeV and rrig = 3, 4, lOGeV are presented at the Table 3. The value of 
af^^y^Mz) is calculated with the error about ±0.004. 
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